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Abstract 

The dynamical instability of a self-gravitating magnetized filamentary cloud was investigated 
while taking account of rotation around its axis. The filamentary cloud of our model is 
supported against self-gravity in part by both a magnetic field and rotation. The density 
distribution in equilibrium was assumed to be a function of the radial distance from the axis, 
Po( r ) = Pc (1 + r 2 /8H 2 )~ 2 , where p c and H are model parameters specifying the density on the 
axis and the length scale, respectively; the magnetic filed was assumed to have both longitudinal 
(z-) and azimuthal (ip-) components with a strength of Bo(r) oc \J po(r). The rotation velocity 
was assumed to be vq^ = f2 c r(l + r 2 /&H 2 )~ l l 2 . We obtained the growth rate and eigenfunction 
numerically for (1) axisymmetric (m = 0) perturbations imposed on a rotating cloud with a 
longitudinal magnetic field, (2) non-axisymmetric (m = 1) perturbations imposed on a rotating 
cloud with a longitudinal magnetic filed, and (3) axisymmetric perturbations imposed on a 
rotating cloud with a helical magnetic field. The fastest growing perturbation is an axisymmetric 
one for all of the model clouds studied. Its wavelength is A max < 11.14 if for a non-rotating cloud 
without a magnetic field, and is shorter when the magnetic filed is stronger and/or the rotation 
is faster. For a rotating cloud without a magnetic filed the most unstable axisymmetric mode 
is excited mainly by self-gravity (the Jeans instability), while the unstable non-axisymmetric 
mode is excited mainly by non- uniform rotation (the Kelvin-Helmholtz instability). The unstable 
non-axisymmetric perturbation corotates with a fluid at r = 2 — 4 H and grows in time. When 
the equilibrium magnetic field is helical, the unstable perturbation grows in time and propagates 
along the axis. A rotating cloud with a helical magnetic field is less unstable than that with a 
longitudinal magnetic field. 
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1. Introduction 

Recent observations have revealed that many molecular clouds contain elongated filamentary 
structures. These filamentary clouds have dense cores as their internal constituents. Most of the 
young stellar objects seem to be associated with dense cores. The evolution, i.e., fragmentation, 
of a filamentary cloud, is therefore interesting in relation to the early phase of star-formation 
processes. 

Filamentary molecular clouds are often associated with magnetic fields which are perpendicular 
to the cloud in the Taurus region (Moneti et al. 1984) and parallel to the cloud in the Ophiuchus 
region (Vrba et al. 1976). It has been suggested that the Orion A cloud has helical magnetic fields 
around the cloud axis (Bally 1989; Uchida et al. 1991). In some filamentary molecular clouds the 
velocity gradient has a component perpendicular to the filament axis, which can be interpreted as 
rotation around the axis (see, e.g., Olano et al. 1988; Uchida et al. 1991; Tatematsu et al. 1993). 
The energies of the magnetic field and the rotation are likely to be comparable to the gravitational 
energy. These magnetic fields and rotation may influence the fragmentation of the filamentary 
clouds. 

Since the early work of Chandrasekhar and Fermi (1953), the linear stability of a cylindrical 
gas cloud has been investigated extensively by Stodolkiewicz (1963), Hansen et al. (1976), 
Nagasawa (1987), Inutsuka and Miyama (1992), and Nakamura et al. (1993, referred to as Paper I 
in the following). These studies have taken account of the effects of magnetic fields, the stiffness 
of the equation of state, the rotation and the collapse of the filamentary cloud in the radial 
direction. However, the stability of a rotating magnetized filamentary cloud against fragmentation 
has not been discussed. Among the studies referred to above, only Hansen et al. (1976) took into 
account the rotation around the axis; however, they considered neither the magnetic fields nor the 
fragmentation in the z-direction. We extend the model of Paper I to include rotation around the 
axis and to discuss the fragmentation of a rotating filamentary cloud with longitudinal or helical 
magnetic fields. 

The model and computation methods are described in section 2. In section 3 a rotating 
filamentary cloud with longitudinal magnetic fields is shown to be unstable against axisymmetric 
and non-axisymmetric perturbations. In section 4 the stability of a rotating filamentary cloud 
with helical magnetic fields is discussed. In section 5 we discuss the application of our stability 
analysis and compare our result to the Kelvin-Helmholtz and Balbus-Hawley instabilities. 

2. Model and Computation Methods 

2.1. Equilibrium Model 

As a model of a filamentary molecular cloud, we considered an infinitely long cylindrical gas 
cloud in equilibrium in which the density and the magnetic field are uniform along the z - axis 
in a cylindrical coordinate system (r, (p, z). The model of Paper I was extended so as to include 
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rotation around the filamentary axis. See, e.g., Bonnell et al. (1992) concerning rotation around 
an arbitary axis. They followed the fragmentation of a rotating finite-length cloud using numerical 
simulations. 

The hydrostatic equilibrium of a filamentary molecular cloud is described by 

dPp d / B 0ip 2 + B 0z 2 ^ | B Qip 2 t „ fVo<p 2 , #(T 
dr dr V 8ir 

and 



-) + ^+»(T + ^)-° w 



Here, P, p, v^, ip, G, B^, and B z are the gas pressure, density, 92-component of the velocity, 
gravitational potential, gravitational constant, and (p- and z-components of the magnetic field, 
respectively. Subscript denotes the quantities in the unperturbed state. The magnetic field of 
the r-component is assumed to vanish in the unperturbed state, So = (0, Bq^, Bq z ). For simplicity, 
the filamentary gas cloud is assumed to be isothermal, 

Po/Po = c s 2 = const. (3) 
A solution satisfying equations (1) through (3) is expressed as 

/ r 2 \-2 

Po = Pc[l + ^pj , (4) 

v = (0,^,0) = {0,rft c (l + r 2 /8i? 2 )- 1/2 ,0} , (5) 

B = (0,B o<p ,B 0z ) , (6) 

Bo* = ^ sin ^( 1 + ^)" 3/2 ' ^ 

Bo, = B c (l + i^Y^f^™^ , (8) 

2 

^0 = 87rGp c iJ 2 ln(l + ^) , (9) 
B 2 

4irGp c H 2 = c 2 + — — (1 + cos 2 #) + 2^1 2 H 2 . (10) 

167T/9C 

All of the symbols with subscript c denote the values at r = 0. The effective radius of this model 
is r = 2V2H. 

The density distribution of our model is the same as those of Stodolkiewicz (1963), Ostriker 
(1964), Nagasawa (1987), and Paper I. Among these models, our solution is the most general in 
the sense that our model incorporates rotation around the axis as well as helical and longitudinal 
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magnetic fields. Our solution has five model parameters: p c , c s , B c , 9, and £l c . The 9 parameter 
denotes the ratio of the ip- and z-components of the magnetic fields and is equal to the pitch 
angle of the magnetic fields at r = oo, 9 = lim tan -1 (Bn^ I Bq z ) . When Q c = 0, our model 
reduces to that of Paper I. When 9 = and f2 c = 0, the magnetic field is parallel to the filament 
axis and the model reduces to that of Stodolkiewicz's (1963) solution for clouds with longitudinal 
magnetic fields. When B c = 0, our solution reduces to Ostriker's (1964) isothermal solution for 
non-magnetized clouds. The sign of B c is taken to be positive in this paper unless otherwise noted. 
In this equilibrium model, the ratios of the magnetic and centrifugal forces to the gas pressure (a 
and P) are spatially constant: 

Bom 2 ~\~ Bq z 2 B<? 
a = — — = j = const (11) 



87rP 8iTp c c s 



and 



= —c " = C onst. (12) 

Cs 2 



The condition for equilibrium [ equation (10)] thus reduces to 

^G Pc H 2 = c s 2 1 + a ( 1 + ™ 2g ) + ■ (13) 

In the following we take H as the unit of length and (2itGp c )~ 1 / 2 as the unit of time. The model 
can thus be specified by three non-dimensional parameters: a, (3, and 9. In the subsequent sections 
we investigate the dependence of the growth rate on a, (3, and 9. 

2.2. Perturbation Equations 

We considered a perturbation of a small amplitude superimposed on the equilibrium cloud 
described above. The linearized equation of motion for the perturbation is given by 

P°\-fr + («o • V)t>i + (wi ■ V)t> ) 

+ VPi + piVV'o + PoV^-i = , (14) 

where v and j = (c/4ir)V x B are the velocity and electric current density, respectively; all of 
the symbols with subscript 1 denote the changes in the quantities due to the perturbation. The 
equation of continuity, the Poisson equation, and the induction equation are respectively expressed 
as follows: 

^ + V • (po»i + Pivo) = , (15) 
A^i = 4ttG Pi (16) 

and 

dB 

= V x (vx x B + v x Bi) . (17) 
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The perturbation is assumed to be isothermal, 



Pi = c s p x . 



(18) 



We obtained the normal-mode solutions of equations (14) through (18), in which all of the physical 
quantities change according to the form 



px(r,<f,z,t) = px(r) exp(— iujt + im<p + ik z z) . 
After some manipulation of equations (14) through (18), it can be rewritten as 



(19) 
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where 



(yi, V2, V3, y4) = (Pi + 



Bi^Bq^ + B\ z Bq z ipoVi r 



PcV'l) Pc9li 



(20) 



(21) 



47T UJ 

The symbol g\ r {= dip\/dr) denotes the change in the gravitational acceleration in the r-direction. 
The elements of matrix A are given explicitly in appendix 1. 

The boundary conditions were set so that the perturbation is regular at r = and 
infinitesimal at r = 00. Since the perturbation equation is singular at r = 0, we obtained regular 
solutions near to r = according to the method of Paper I. The asymptotic solutions in the 
region r » AH were also obtained according to Paper I. 

2.3. Numerical Procedure 

The eigenvalues and eigenfunctions were obtained essentially with the methods of Paper 
I. Namely, we integrated the perturbation equation both from r = and 00 for a given 
lo (= Lo r + iuji) and checked whether the integrated solutions could be connected continuously at 
a middle point. When a solution satisfied the boundary conditions at r = and 00 for a given u, 
we regarded the given lo as beeing an eigenvalue. 

We searched for eigenvalues using the bisection method. When f2 = 0orm = # = 0, the 
eigenvalue (u> 2 ) is always real and can be obtained by the usual bisection method. When otherwise, 
we extended the bisection method for a complex eigenvalue. The technical details concerning the 
extended bisection method are given in appendix 2. 

Equation (20) has a singularity at vq^ = ruj r /m when m / 0. The result depends on the 
integration path (see, e.g., Lin 1945a, b; Kato 1987). We integrated equation (20) over r along the 
real axis. This integration path is correct only for unstable modes. Although this integration path 
is not correct for damped modes, we applied it for simplicity, since we are not much interested in 
damped modes. 
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3. Rotating Clouds with Longitudinal Magnetic Fields 

In this section we consider the stability of a rotating cloud with longitudinal magnetic fields. 
The parameter 9 is fixed to be 0° in this section, except where otherwise stated. 

3.1. Sausage (m = 0) Mode 

In the beginning we consider the sausage mode (m = 0) instability of a rotating non- 
magnetized filamentary cloud. The eigenvalue (w) is purely imaginary for the unstable mode. 
Figure 1 shows the dispersion relation for the model of (a, (3) = (0, 1) by the thick curve. 
The abscissa and ordinate are the non-dimensional wavenumber (k z H) and the non-dimensional 
growth rate (uJi/^2iTGp c ), respectively. There is only one unstable sausage mode for the models 
with o = 0, although there are some unstable sausage modes for a > (Paper I). The sausage 
mode is unstable only when the wavenumber is smaller than a critical one, k Z)CT H = 0.960. The 
growth rate has its maximum, a^ imax = 0.606 \J2ixGp c ^ at k z H = 0.467. The dashed and thin 
curves denote the growth rate of the most unstable mode for (a, [3) = (0, 0) and (a, [3) = (1, 0), 
respectively. The rotation as well as the longitudinal magnetic fields increase the maximum 
growth rate (= o^max), the wavenumber of the fastest growing perturbation (= fc 2;max ), and the 
critical wavenumber (= k Z)Cr ). This is because both rotation around the axis and longitudinal 
magnetic fields support the gas against gravity in the radial direction, but does not operate in 
the z-direction. The gas temperature is lower for fixed p c and H when a and (3 are larger [ see 
equation (10) ]. Accordingly the Jeans length, the typical length scale for fragmentation due to 
self-gravity, is shorter. The wavenumber of the fastest growing perturbation is almost the same for 
(a, (3) = (0, 1) and (1, 0). This implies that k Z:inax is a function of c s for = 0°. The maximum 
growth rate is smaller for (a, (3) = (0, 1) than for (a, [3) = (1, 0). This means that the Parker 
instability increases the growth rate when a magnetic field is present. 

Figure 2 shows the dependence of the growth rate on (3 for a = 0. As (3 increases, w^max) 
fc 2]mEK , and k z ^ cr increase. The wavenumber of the the fastest growing perturbation can be 
approximated by 

k z , max H = 0.72 [(1 + a + (3) 1 ^ - 0.6] (22) 

for 6 = 0°. Figure 3 shows the accuracy of equation (22). The upper panel shows k Z)inax as a 
function of a + (3. Equation (22) is drawn by the curve and the numerically obtained values 
(fez, max) are plotted with filled circles. The lower panel shows the deviation from equation (22). 
Equation (22) gives a good estimate for k z ^ mSuX with an error of 2% for < a + (3 < 6. Equation 
(22) is used for a comparison with a filamentary molecular cloud, Orion A, by Hanawa et al. 
(1993). See Hanawa et al. (1993) for an application of equation (22). The approximate dispersion 
relation is given by 



2 , n k z H 0.89 + 1.4a 2 2 

: " - 4nGpc TTkJI 1 + 1.25a + Cs h 

k 7 H 0.89 + 1.4a H 2 k 7 2 




1 + k z H 1 + 1.25a 1 + a + (3 



(23) 
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for the most unstable sausage mode for all the models of 6 = 0° . The upper and lower expressions 
of equation (23) are equivalent [see equation (13)]. The second term on the right-hand side of 
equation (23) denotes the stabilization due to thermal pressure, i.e., the dispersion relation of the 
sound wave. When a and are larger, c s is lower, and, accordingly, the growth rate is larger. The 
maximum growth rate obtained from equation (23) agrees with the numerical results with an error 
of 5% for 0.5 < a < 10 and < < 10 , and for a = and < < 2. 

Figure 4 shows the fastest growing perturbation for (a, 0) = (0, 1). The density distribution 
in the r — z plane is indicated by the grey scale and the velocity field is denoted by arrows. The 
amplitude of the perturbation is taken to be pi/po = £ cos ik z z) at r = 0. To emphasize the 
density contrast we have taken e = 0.65 in figure 4. The relative density perturbation (|pi/po|) 
has its maximum on the axis (r = 0). The velocity perturbation is dominant in the z-direction. 
The velocity perturbation in the z-direction has its maximum, v± z = — 0.941 ec s sin (k z z), on the 
axis and monotonically decreases as r increases. The velocity perturbation in the r-direction has 
its maximum, v\ r = — 0.128 ec s cos (k z z), at r = 2.73 H. The 92-component of the velocity is 
shown in the right panel. The solid and dashed curves denote the velocity in the perturbed state 
and in equilibrium, respectively. The change in the rotation velocity is small with a maximum, 
viip = 0.193 ec s cos (k z z), at r = 1.76 H. 

The eigenfunction shown in figure 4 is similar to that of the the fastest growing perturbation 
for a magnetized non-rotating cloud (see figure 8 of Paper I). A filamentary cloud supported in 
part by either rotation or longitudinal magnetic fields fragments mainly in the z-direction. We 
expected in Paper I that a disk perpendicular to the rotation axis is formed by the nonlinear 
growth of the perturbation. Tomisaka (1993) and Nakamura et al. (1993, private communication) 
showed this by numerical simulations for magnetized filamentary clouds. We also followed the 
nonlinear evolution of a rotating filamentary cloud with numerical simulations, and confirmed the 
formation of a small disk perpendicular to the axis. The results of these numerical simulations 
will be reported in a future paper. 

Figure 5 is the same as figure 4, except for (a, 0) = (0.5, 0.5). The velocity perturbation 
is dominant in the z-direction. The velocity perturbation in the z-direction has its maximum, 
viz = — 0.886ec s sin(/c 2 z), on the axis when p\/po = £ cos(k z z). The velocity perturbation 
in the r-direction has its maximum, v\ r = — 0.175ec s cos(k z z), at r = 2.87 H. The change in 
the rotation velocity is at most vi^ = 0.101 ec s cos (k z z) at r = 1.63 H. The characteristics of 
the most unstable sausage mode for (a, 0) = (0.5, 0.5) are between those for (a, 0) = (1, 0) 
and (0, 1), except for the generation of B^. When a0 7^ 0, the magnetic field is twisted due 
to a change in the rotation velocity [B\ v = — 0.249 e sin (k z z) at r = 1.51 H, at most while 
B 0z = 3.545(1 + ^/m 2 )' 1 ]. 

3.2. Kink (m = 1) Mode 

Rotating filamentary clouds are also unstable against non-axisymmetric perturbations. We 
found five modes, three of which are pure oscillations (a;, = 0) for any k z . The remaining two 
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modes have complex conjugate eigenfrequencies. Figure 6 shows the dispersion relation of the 
unstable kink mode (m = 1) for a = and (3 = 1. The real and imaginary parts of the 
eigenfrequency are denoted by the dashed and solid curves, respectively. The dispersion relation 
is symmetric with respect to k z , i.e., uj(— k z ) = oj{k z ). We thus restrict ourselves to k z > in 
the following. The three horizontal dotted lines in the figure denote the angular velocities , S7, 
at r = 2, 2\/2, and 4, respectively. The unstable modes have a corotation point where the phase 
velocity coincides with the fluid velocity [oj r /m = fi(r)]. The corotation point is located near to 
the effective radius, r = 2y/2H, when the growth rate is large. The growth rate is very small 
{uji < 1 x 1(T 4 y/2-KGpc) in the region \Hk z \ < 0.04, and large in the region 0.1 ^ Hk z < 0.477. 
Because of the analytical and numerical difficulty at the corotation point [uj r = mfi(r)], we could 
not find whether u>i = in the region \Hk z \ < 0.04. When k z is larger than the critical value 
(& 2 ,cr = 0.477 H~ l ), the growth rate is again small (uji < 10~ 2 ^2irGp c ). 

The mode having m = —1 is physically the same as that having m = 1, except for the sign 
of Lo r . When m = —1, a mode with uj t < becomes unstable and a mode with to r > is a pure 
oscillation. Only when the fluid velocity (= fi) and the phase velocity of the wave (= uj r /m) have 
the same sign can the mode be unstable. 

Figures 7 and 8 are the same as figure 6, except for ( a, (3) = (0, 0.5) and (0, 2), respectively. 
The growth rate is larger when (3 is larger. The non-rotating filamentary cloud without magnetic 
fields (a = (3 = 0) has two neutrally stable kink modes with positive and negative real frequencies 
(see also Nagasawa 1987). Only the mode having a corotation point becomes unstable; the other 
mode remains pure oscillatory. We can thus conclude that this instability is due to the resonance 
of the wave to the fluid rotation, and, accordingly, is of the Kelvin-Helmholtz type. Note that 
non-rotating magnetized filamentary clouds are unstable against the kink mode because of the 
Parker instability (Paper I). 

When \k z \ is small, the growth rate is very small (even if it is positive) and the oscillation 
frequency is almost independent of (3. For long-wavelength perturbations v z has a large amplitude 
in the region very far from the axis (r ;> 10 H). For \k z \ > k Z)Cr the growth rate of the kink 
mode is very small (u>i < 1 x 10 -2 ). The critical wavenumber is Hk z>CT = 0.301, 0.477, and 
0.788 for (a, [3) = (0, 0.5), (0, 1), and (0, 2), respectively. The growth rate has its maximum 
[cj i / v /27rG>c = 0.0564, 0.135, and 0.255] at Hk z = 0.225, 0.323, and 0.478, for (a, f3) = (0, 0.5), 
(0, 1), and (0, 2), respectively. 

Figure 9 is the same as figure 4, except for the most unstable kink perturbation for 
(a, (3) = (0, 1). Figures 9(a) and (b) denote the cross section in the r — z plane and that of 
z = 0, respectively. The growth rate and wave number are uji = 0.135 \/2irGp c and Hk z = 0.323, 
respectively. The amplitude of the perturbation is normalized so that the density perturbation is 
at most pi/ po = 0.5. When the kink mode is excited, the density ridge is twisted, as can be seen 
in figure 9. 

Figure 10 is the same as figure 6, except for (a, (3) = (0.5, 0.5). This model is intermediate 
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between those of (a, (3, 9) = (1, 0, 0°) and (0, 1, 0°), the former of which is excited by the Parker 
instability. For a fixed value of a + (3 the kink mode has a larger growth rate when a is larger. 

4. Rotating Cloud with Helical Magnetic Fields 

In this section we discuss the sausage mode instability of rotating filamentary clouds with 
helical magnetic fields (a/36 / 0). 

Figure 11 shows the dispersion relation of the most unstable sausage mode for 9 = 30°, 60°, 
and 90°, while a = (3 = 0.5. The real (figure 11a) and imaginary (figure lib) parts of the 
eigenfrequency are shown as a function of k z . The sausage mode grows and propagates with a 
phase speed of to r /k z in the z-direction. The eigenfrequency of the damped mode is the complex 
conjugate of that of the unstable mode. The eigenfrequency is antisymmetric with respect to k z , 
i.e., u)(—k z ) = —u)(k z ). The phase velocity is thus always positive for both the unstable and 
damped modes. When 9 < 0, the real part of the eigenfrequency is negative, uj r < 0. Note that 
O c is taken to be positive in our computations. The phase velocity is to r /k z ~ 0.025 H\J2irGp c 
for 9 = 30°. The phase velocity is faster for a larger 9. 

As 9 increases, w !]max , k z ^ m3iX , and /c ZjCr decrease, and, accordingly, the cloud becomes less 
unstable. This is due to an increases in the sound speed, c s , for fixed a and (3 [ see equation (13)]. 
Both when the magnetic field is longitudinal and helical, the most unstable mode is mainly excited 
by a self-gravitational instability. The eigenfrequency oj is complex for [39 / 0, while it is pure 
imaginary or real for f39 = 0. The perturbation therefore grows in time and propagates in the 
z-direction only when the cloud rotates and is threaded by a helical magnetic field. 

Figure 12 shows the fastest growing perturbation for (a, (3, 9) = (0.5, 0.5,60°). Figure 12 
is similar to figure 5, except for the phase shift in the z-direction. The density perturbation is 
proportional to p\ oc cos(k z z — uj r t + 5 P ), where the phase shift, 5 P , is a function of r, and is 
taken to be 5 P = at r = 0. The phase shift is large for the azimuthal components of the 
velocity and magnetic fields, e.g., v\ v = 0.141 e cos (k z z — uj r t + 0.561) at r = 1.73 H and 
Bip = 0.682esm(k z z — u) r t + 0.356) at r = 1.25 H. The upper panel of figure 12 shows the 
z-dependence of vi^ (the dashed curve) at r = 1.73 H and Bi^ (the solid curve) at r = 1.25 H, 
both of which have their maximum amplitudes there. 

The phase shifts are related to the propagation of the wave. Consider a velocity perturbation 
proportional to v\ r oc cos (k z z) and v\ z oc sin (k z z). The velocity perturbation changes v v and 

according to the azimuthal components of the equation of motion and the induction equation. 
When 9 = 0, the changes in the azimuthal components have the dependence B\^ oc sin (k z z) and 

oc cos(k z z) (see section 3). When (3 = 0, they have the dependence B\^ oc cos(k z z) and 
viip oc sin (k z z) (see Paper I). When (39 / 0, both types of the above-mentioned perturbations 
are produced. These changes in v v and B v modify v r and v z according to the equation of motion. 
As a result, velocity perturbations given by v\ r oc •&m.(k z z) and v\ z oc sm(k z z) are also induced 
and the velocity perturbations apparently propagate. 
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5. Discussion 

As shown in the previous sections, a rotating magnetized filamentary cloud suffers from 
various instabilities. This is in part because our model includes rotation as well as magnetic fields 
as well as the self-gravity of the gas. In this section we discuss the relationship between our 
model and other theoretical studies conserning the instability of a rotating gas cloud. In the last 
paragraph we also discuss the application of our model to TMC 1 (Taurus Molecular Cloud 1). 

5.1. Relationship to the Kelvin-Helmholtz Instability and the Balbus-Hawley Mechanism 

Since the angular velocity is non-uniform in our model, we suppose that the Kelvin-Helmholtz 
instability may be involved in the stability of our model cloud. As discussed by Glatzel (1987a,b), 
the resonant interaction of two neutral modes produces a couple of growing and damped modes 
and causes an angular momentum transfer across the corotation point. These characteristics are 
shared with the unstable kink mode for a rotating cloud without a magnetic field. As shown in 
subsection 3.2, the growing and damped modes appear in pairs. We confirmed that the angular 
momentum is transferred from the region inside the corotation point to that outside the corotation 
point when the unstable kink mode is excited. In the growing mode, p\ and v\ v are anti-correlated 
(pi* ■ viu, < 0) inside the corotation point and are positively correlated (p\* ■ v\ v ) outside the 
corotation point, where pi* denotes the complex conjugate of p\. This means that the central 
region loses angular momentum and the outer part receives it. The angular-momentum transfer 
produces an energy excess, which causes the perturbation to grow. 

The energy excess due to the angular momentum transfer is proportional to the difference 
in the angular velocity across the corotation point. Correspondingly, the growth rate of the 
Kelvin-Helmholtz instability is bounded by the angular velocity gradient, uj; l < max \(\/2)rd£l/dr\ 
(Sung 1974; Hanawa 1986, 1987; Fujimoto 1987), when the self-gravity of the gas and magnetic 
fields are negligibly small. In our model, the angular velocity gradient, 



has its maximum, 0.385 Sl c at r = 4H, and a somewhat smaller value 0.355 Q c at the effective 
radius, r = 2y/2H. It may be not a chance coincidence that the growth rate of the kink mode is 
large when the angular velocity gradient is large at the corotation point. 

The unstable kink mode might be excited in part by another effect related to a large (3. 
When (3 is large, the sound speed (c s ) and thermal pressure are low [ see equation (13)]. The 
low gas pressure strengthens the effect of the self-gravity relatively and suppresses the Jeans 
instability less effectively. This effect is appreciable in the axisymmetric mode when (3 ^ 1. To 
investigate this effect on the kink mode, we made an experimental model in which the right-hand 
side of equation (16) was multiplied by a factor of < e < 1, A^i = 4nGp\ . By decreasing e 
from unity, the change in the gravitational potential was reduced artificially. As e decreases, u> r 
increases and u>i decreases. The increase in u> r is ascribed to a reduction in the self-gravity, since 




(24) 
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the self-gravity lowers the propagation of the sound wave. The increase in uj r shifts the corotation 
point inwards, and, accordingly, decreases the angular velocity gradient thereof. The decrease in 
oji may be in part due to a decrease in the angular velocity gradient at the corotation point. Thus, 
unfortunately, we could not decouple the effect of low thermal pressure on the Jeans instability 
from other effects of large f3. 

When a[3 / 0, the unstable sausage mode is in part due to the Balbus-Hawley (1991) 
mechanism, although it is mainly due to the Jeans instability. Balbus and Hawley showed that a 
differentially rotating disk is highly unstable if it is threaded by a vertical magnetic filed. In such 
a disk, an outwardly displaced fluid element tries to enforce a rigid rotation, and thus rotates too 
fast for its new radial location. Then, the excess centrifugal force drives the element still faster 
outward. As can be seen in subsection 3.1, the unstable sausage mode generates B r and B^. The 
twisted magnetic fields transfer angular momentum along the field line. We confirmed that the 
fluid element displaced inward loses angular momentum due to magnetic drag when the unstable 
sausage mode is excited for a(3 / 0. Because of the angular momentum gain, the fluid element 
contracts further due to a decrease in the centrifugal force. In order to evaluate the contribution 
of the Balbus-Hawley mechanism to the unstable sausage mode, we again computed the artificial 
models in which the change in the gravitational potential was attenuated by a factor of e. The 
growth rate, LOi, decreases along with a decrease in e, and diminishes at s = 0. We thus conclude 
that the unstable sausage mode is excited mainly by self-gravity. When e = 0, the increase in the 
magnetic pressure near to the axis (r = 0) stabilizes our model cloud against the Balbus-Hawley 
mechanism for m = 0. We confirmed that the sausage mode is stabilized against the Parker 
instability by the same mechanism for (a, (3) = (1, 0) when e = 0. 

5.2. Comparison with Habe et al.'s Simulation 

Here, we comment on the collapse of a self-gravitating cloud triggered by torsional Alfven 
Waves. Habe et al. (1991) showed by a numerical simulation that a cloud with helical magnetic 
fields collapse more easily than one with longitudinal magnetic fields. In their simulation, helical 
magnetic fields are generated by rotation and cause a rotating cloud to collapse due to magnetic 
pinching. It seems as if our stability analysis contradicts their simulation, since a rotating cloud 
with helical magnetic fields is less unstable than that with longitudinal magnetic fields in our 
analysis. The apparent contradiction comes from a difference in the initial models. While our 
initial model is in equilibrium, Habe et al's (1991) started their numerical simulation from a cloud 
threaded by helical magnetic fields. Their initial model is not in equilibrium. And our stability 
analysis cannot be applied directly to their initial model. 

5.3. Application to TMC 1 

Finally, we discuss the application of our model to a filamentary molecular cloud. TMC 1 
contains 5 dense cores which are more or less regularly spaced with an average separation of 4'. 
Since the apparent diameter is 2', the ratio of the average separation to the filament diameter 
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is estimated to be 2. Substituting this value into our model, we obtain Hk max ~ 0.63. This 
means that the magnetic and/or centrifugal forces are comparable to the gravitational force 
(a + f3 ~ 1). This is consistent with the velocity structure of TMC 1. Olano et al. (1988) found 
a velocity gradient across the filament axis, which can be interpreted as rotation around the axis. 
The rotation velocity is comparable to the velocity dispersion. This implies that the cloud is 
supported against gravity in part by rotation, and that the centrifugal force is as strong as the 
turbulent pressure, ((3 « 1). 
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comments. This work is financially supported in part by the Grant-in- Aid for General Scientific 
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Appendix 1. The Matrix Elements of A 



Matrix A is expressed as 
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Appendix 2. The Bisection Method for Complex Eigenvalues 

In this appendix we describe the method used to search for a complex eigen value. According 
to Paper I, the condition for an eigenvalue is expressed as 



X(w) = det 
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\yl\r;uj) yf\r; uj) yf\r; uj) y^\r; uj) ) 



where [yi \ y%\ yi\ yf 1 ] for i = 1, 2, 3, and 4 denote independent solutions satisfying either of 
the boundary conditions at r = or 00. We first evaluate x( w ) f° r ^ = w i> ^2, and ^3- The 
points uj = o>2j an d w 3 are on the vertexes of a right isosceles triangle on the complex uj plane 
[ see figure (13)]. We next compute 



77(^1 w 2 w 3 ) = log[x(^i)/x(w 2 )] + log[x(w 2 )/x(w 3 )] + log[x(w 3 )/x(^i)] , 
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where each logarithm operation takes the principal value. The value of t](uji u>2 UJ3) can take 
either or ±2iri. When 77 = 2m, there is a zero point of x( w ) inside the triangle. Note that 

V = H 1 /x)dx- 

When 77 = 2-7TZ for the triangle u>\ — u>2 — u>s, we integrate the perturbation equation for 
UJ4 = (u>2 + ujs)/2 and compute x( w 4)- Either rj(ui, w 2 , L04), rj(uji, uj^, us), or 77(0^3, W4, w 2 ) is 
identical to r](uJi, cj 2 , W3), 



7?((Jl, W 2 , W 3 ) = 7?(o;i, U> 2 , UJ4) + 7?(^l, W4, UJ 3 ) + 77(^3, W 4 , W 2 ) . 



(76) 



In most cases, we obtain 77(0^3, lj 4 , w 2 ) = 0. [ The frequency of 77(^3, w 4 , w 2 ) / is very small 
when the triangle uj\ — w 2 — W3 is small ]. Then, the region containing the zero point of x (v) = is 
reduced to half in area. The reduced region is again a right isosceles triangle and can be squeezed 
by successive reduction. Finally, we can obtain an eigenvalue, uj, with a sufficiently small error. 
Each iteration makes the estimated error smaller by a factor of y/2 in this extended bisection 
method. 
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Fig. 1. — Dispersion relation for models with (a, (3) = (0, 1). The thick curve denotes the growth 
rate, u>i, as a function of the wavenumber, k z H. The thin and dashed curves denote the growth rate 
of the most unstable mode with (a, (3, 9) = (1, 0, 0°) and (0, 0, 0°), respectively, for comparison. 

Fig. 2. — Dependence of the growth rate, Wj, of the sausage (m = 0) mode on (3 for a = 0. Each 
curve denotes the growth rate of the unstable sausage mode with (a, (3) = (0, 0), (0, 1), and (0, 2). 
The ordinate is the growth rate, Ui, and the wave number, k z . 

Fig. 3. — Wave number of the most unstable perturbation as a function of a + (3 for 9 = 0°. The 
ordinate and abscissa of the upper panel are the wavenumbers of the most unstable perturbation, 
kmaxH and a + (3, respectively. The filled circles denote the numerically obtained data points, and 
the curve denotes the fitting formula, equation (22). The lower panel shows the error of the fitting 
formula. The ordinate is (error) = (fitting formula) /(numerical data) — 1. 

Fig. 4. — Cross section of the model filamentary cloud perturbed by the fastest growing sausage 
mode (main panel). The model parameters of the equilibrium model are (a, (3, 9) = (0, 1, 0°). The 
abscissa is the z-axis and the ordinate is the radial direction. The density is indicated by greyness, 
the scale of which is shown on the left side of the panel. The arrows denote the velocity field on 
the r — z plane. The right-hand panel shows Vu> as a function of r at z = 0.0. The dashed and 
solid curves denote the values in equilibrium and in the perturbed state, respectively. The growth 
rate and the wave number of this perturbation are a>j = 0.606 and k z = 0.467, respectively. 

Fig. 5. — Same as figure 4, but for (a, (3, 9) = (0.5, 0.5, 0°). The growth rate and the wave number 
of this perturbation are u>i = 0.640 and k z = 0.478, respectively. 

Fig. 6. — Dispersion relation of the kink (m = 1) mode for (a, (3) = (0, 1). The dashed and 
solid curves denote u r and Wj, the oscillation frequency and growth rate, respectively. The three 
horizontal dotted lines denote the angular velocities, Q, at r = 2, 2^/2, and 4, respectively. 

Fig. 7. — Same as figure 6, except for (a, (3) = (0, 0.5). 

Fig. 8. — Same as figure 6, except for (a, (3) = (0, 2). 

Fig. 9. — Same as figure 4, except for the most unstable kink mode with (a, (3) = (0, 1). The 
wave number and growth rate are Hk z = 0.323 and Wj jmax = 0.135 ^/(2irGp c ), respectively. The 
amplitude of the perturbation is normalized so that pi/p c = 0.5 at maximum, (a) The cross 
section in the r — z plane, (b) The cross section of z = . 

Fig. 10. — Same as figure 6, except for the kink mode of (a, (3, 9) = (0.5, 0.5, 0°). 
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Fig. 11. — Dispersion relation of the sausage mode with (a, (3) = (0.5, 0.5). The curves denote 
the growth rate as a function of k z for 9 = 30°, 60°, and 90°. Figures 11 (a) and (b) show uj r and 
LOi, respectively. The dotted curve denotes uji for a model with 9 = 0° for comparison. 



Fig. 12. — Same as figure 4, except for the sausage mode with (a, (3, 9) = (0.5, 0.5, 60°). The 
growth rate and wave number are a; = (0.0199 + 0.594 i) ^/(27rGp c ) and k z H = 0.432, respectively. 



Fig. 13. — Search for a complex eigenvalue on the complex u plane. The region for the search 
narrows by a factor of 2 in area by each iteration. 



